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Compartmentalization into biochemically distinct organelles constantly exchanging material is 
one of the hallmarks of eukaryotic cells. In the most naive picture of inter-organelle transport driven 
by concentration gradients, concentration differences between organelles should relax. We determine 
the conditions under which cooperative transport, i.e. based on molecular recognition, allows 
for the existence and maintenance of distinct organelle identities. Cooperative transport is also 
shown to control the flux of material transiting through a compartmentalized system, dramatically 
increasing the transit time under high incoming flux. By including chemical processing of the 
transported species, we show that this property provides a strong functional advantage to a system 
responsible for protein maturation and sorting. 



Eukaryotic cells contain many specialized compart- 
ments (organelles) constantly exchanging molecules 
through complex energy-consuming transport processes. 
Along the secretory pathway for instance, proteins 
and lipids synthesized in the endoplasmic reticulum 
(ER) are transported to the Golgi complex (itself sub- 
compartmentalized into stacked cisternae) for matura- 
tion and sorting, before being conveyed to particular cel- 
lular locations pQ. ER and Golgi (and even the different 
Golgi cisternae) are characterized by fairly different lipid 
and protein compositions [2]. There is a widespread in- 
terest in understanding the molecular and physical bases 
permitting these organelles to maintain their identity, 
namely their specific chemical composition, while con- 
stantly exchanging material [3]. 

In a system where fluxes are linearly related to concen- 
tration differences {i.e. satisfying Fickian diffusion), sta- 
tionary concentration gradients can only be maintained 
by external fluxes. Fluctuations of the fluxes yield fluc- 
tuations of the local composition, and robust compart- 
ment identity (namely the existence of stable concentra- 
tion heterogeneity) is not to be expected. In the secre- 
tory pathway however, transport is heavily influenced 
by molecular recognition: transported molecules, carri- 
ers and recipient organelles cooperate through complex 
networks of molecular interactions [H |S). The goal of 
this paper is to study the consequences of such cooper- 
ative transport on the stationary state of a simple (2- 
compartment) system, and on its possible function in 
protein maturation and sorting. Our model generalizes 
and extends the vesicular transport model proposed by 
Heinrich and Rapoport [6], where the rates of vesicle 
exchange between compartments are influenced by their 
composition. We discuss below both the case of a closed 
system and of an open system with incoming and outgo- 
ing fluxes. By including chemical transformation of the 
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transported species, we show that cooperative transport 
can strongly increase the accuracy of a system responsi- 
ble for protein maturation and sorting. 



I. STATIONARY COMPARTMENT 
DIFFERENTIATION IN A CLOSED SYSTEM 

We first consider a single protein species distributed 
between two compartments constantly exchanging mate- 
rial, and indicate how these results might be extended to 
a multicomponent system. 



A. One species system 

We assume below that the total mass of the system 
(and the mass of each compartment) is maintained con- 
stant by an unspecified regulatory mechanism, so that 
the evolution of the concentration C\ in compartment 1 
can generically be described by the Master equation [TJ: 

dtd = h- J^ 2 (Ci,C 2 ) + J 2 ^i(C 2 ,d) (1) 

with a similar expression for compartment 2 obtained 
by the transformation 1 ■<-» 2. Here, J a ^,p(C ai Cp) is 
the mean flux from compartment a (with concentration 
C a ) to compartment j3 (of concentrations Cp). Compart- 
ments will naturally reach different concentrations if they 
follow different exchange rules. We focus on the more in- 
teresting case where J\-^ 2 {C\, C 2 ) = J2-n(Ci, C2). The 
source and sink term l\ in Eq[T] may include both ex- 
ternal fluxes in and out of compartment 1 and chemical 
transformation within this compartment. 

The transport of cargo between organelles may be sep- 
arated into three distinct steps; step 1: cargo packaging 
inside a membrane-based carrier, such as a small protein- 
coated vesicle or a membrane tubule [8 , step 2: the ac- 
tual transport between secreting and receiving compart- 
ments, often involving molecular motors moving along 
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FIG. 1. (a) Location of the critical region in the parameter space {C s /C to t, C//C to t} where stationary compartment differen- 
tiation occurs in a closed system (shaded blue, Eqj5|. Increasing the total concentration Ctot moves the system along the red 
line (arrow), (b): Variation of the stationary compositions (black) and flux (red) with the total concentration, showing the 
breaking of symmetry for Ctot > C* ot (Eq[5} . 



cytoskeletal filaments [51 [TP], and step 3: fusion of the car- 
rier with the receiving organelle (see [TT] for a review). 
For step 1, we call J s the total flux of material secreted by 
a compartment, and S the fraction of this flux (a number 
between and 1) occupied by the species of interest. Let 
us assume step 2 to be infinitely fast (we relax this hy- 
pothesis in Sec |I C[ ). In this case, any vesicle secreted by 
compartment 1 immediately merges (step 3) either with 
compartment 2, with a probability Pi->2, or back with 
compartment 1 (with probability Pi->i = 1— Pi_>2)- The 
mean flux from compartment 1 to compartment 2 may 
thus be written: 

Jx^2{Cx,C 2 ) = J s {C 1 )S(C 1 )P l ^{C u C 2 ). (2) 

For a closed system with fixed concentration C to t — 
C\ + Ci (no source and sink term), the symmetric 
state: C\ — C2 — C to t/2 is a stationary solution 
(dtCi = d t C2 = 0). Linear stability analysis [7] 
shows that the symmetric solution is unstable provided 
=c 2 =c tot /2 ^ 0' The case of particles ran 
domly entering transport vesicles which are secreted at 
constant rate and fuse non-specifically with either com- 
partment corresponds to a linear flux (</i->2 ~ Ci) (akin 
to passive diffusion) and leads to identical compartments. 
Spontaneous compartment differentiation can only occur 
if the flux reaching the second compartment decreases 
with increasing concentration in the first one, and this 
requires non- linear transport (cooperativity) . 

Let us assume the fluxes to have two rather univer- 
sal types of non-linearity as a function of concentration. 
Firstly, the out-going flux of a given species should sat- 
urate at high species concentrations. This can be due 
to the limited capacity of transport vesicles, the limited 
availability of vesicle-coating proteins, or the formation 
of aggregates inapt for transport in a compartment be- 
yond a critical concentration. For simplicity, we choose 
to keep the flux of secreted vesicles constant (and write 
it J s = K C S ), although direct interactions between car- 
goes and coat proteins are known to exist [5]. The pack- 



aging fraction S is assumed to saturate beyond a concen- 
tration C s following a Michaelis-Menten saturation [T2"] : 



S(Ci 



Ci + C s 



(3) 



Secondly, vesicle fusion is known to be strongly reg- 
ulated by specific molecular interactions, including, but 
not restricted to, interactions between matching pairs of 
SNAREs[4 . Quantitative models have shown the impor- 
tance of this step for the generation and maintenance of 
non-identical compartments, using fairly detailed mathe- 
matical modeling of the different pairs of SNAREs ]H |T3] 
and/or extensive numerical simulations [14J . Numerous 
factors can however influence the delivery of transport 
vesicles, including specific interactions between the cargo 
and molecular motors [5]. Here, we adopt a very generic 
treatment of specific fusion, where the fusion probability 
P\^2 deviates from its nonspecific value because of two- 
body interactions between constituents of the vesicle and 
the receiving compartment: P±^2 — 1/2 ~ S(Ci)C2- Af- 
ter normalization, the probability may be written: 



Pi- 



C f + S(Ci)C 2 



2C / + S(C 1 )(C 1 + C 2 ) 



(4) 



where Cf is the typical concentration beyond which 
specific fusion becomes relevant. Within the descrip- 
tion outlined in Eqs.(3p), linear transport corresponds 
to both characteristic concentrations being very large: 
C s ,Cf ^> Ctot • 

Spontaneous symmetry breaking (enrichment of one 
compartment at the expense of the other) occurs when 



■2>C tot /2 



< 0. As shown in Fig la, this always 



happens at high enough concentration 
C* t J = AC s C f {C s + Ct 



> C* ot , with 



(•5) 



Beyond this threshold, any small perturbation from the 
symmetric state brings the compartments into a stable 
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asymmetric steady-state. As a consequence, the concen- 
tration of the least concentrated compartment (compart- 
ment 2, say) and the flux Ji_>.2 of material exchanged be- 
tween compartments both decrease with increasing con- 
centration when Ctot > Cj ot , Fig[l]}. At high concentra- 
tion, the asymptotic solution reads C2 ~ 2C/C s /Ct t- 

Although the actual location of the critical line de- 
fined by Eq|5] depends on the model (Eqs.(3p|) for the 
exchange flux Ji_>2 (Eq[2]), its existence does not. This 
critical behaviour is very general and stems from the pres- 
ence of two competing effects : cooperative fusion pro- 
motes protein enrichment (and increases with decreas- 
ing Cf), while saturation of protein packaging (beyond 
a composition C s ) limits transport. Including the pres- 
ence of different types of coat and fusion proteins does 
not fundamentally alter this picture [6]. 



B. Extension to a n-species system 

Extending the analysis presented above to a n- 
component system is rather straightforward. Let us call 
C l a the concentration of the species i in the compartment 
a (a = 1, 2). The concentration of all species in compart- 
ment a can be defined as a vector C a = [C*, (7^, C™], 
and satisfies the Master equation: 



d t ci 



p - p 



J, 



(6) 



where J l a ^,p is the mean flux of the species i from the 
compartment a to the compartment /3, and P a is a net 
source and sink term including both the presence of ex- 
ternal fluxes of species i in and out of compartment a, 
and chemical transformation involving species i in com- 
partment a. 

For a closed system (no source and sink term), the total 
concentration for the i-th specie is fixed: Cl ot — C l a + Cy. 
All the equations may thus be written for the fractions 
C = Q/Qot, satisfying c/>\ = 1. Then <f> 2 = 
becomes implicit and the master equation is now written 
only as a function of <p = <fii '■ 

dtP = -jU 2 (<l>, 1 - <t>) - 0) (7) 



with the normalized fluxes ?J 

■J Ct— 



Assum- 



ing as before that both compartments follow identical 
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is a sta- 



exchange rules, 4>i/2 = 1 ~ 01/2 

tionary solution. The linear stability of the symmetric 
solution is determined by the Jacobian matrix M : 



Mi. 



-2 {d^t 



'4>i 



/2 



(8) 



The symmetric state is unstable, and spontaneously 
evolves towards a non-symmetric state if If M has at 
least one positive eigenvalue. 

In a multi-component system, the fluxes can be written 
similarly to the main text: 



The functions J a , S l a and P a ^p may contain various non- 
linearities. In particular P a ^/3 may involve any combi- 
nation of pair interactions {S l a ,C^} which can lead to 
a very rich behaviour. One could in particular describe 
in this way the transport of proteins directly interacting 
with the secretion (coat proteins) or the fusion (SNAREs) 
machinery. 



C. Influence of a finite vesicle fusion time 

If vesicular transport between secreting and receiv- 
ing compartment (the so-called step 2 in I A I is not in- 
finitely fast, vesicles will dwell for some time the inter- 
compartment region, and will have a complex distri- 
bution of concentration, reflecting the concentration of 
the emitting compartment at the time of their secre- 
tion. While this situation appears much more complex 
than the one described above, we now show, restricting 
ourselves to a one-species system for simplicity, that a 
model with inter-compartment dwelling of vesicles can 
be mapped to the simpler model with immediate fusion 
of vesicles. 

Each vesicle can carry a given amount of proteins, and 
a vesicle budding from or merging with a compartment 
will change the concentration of this compartment. Let 
us call C v the resulting change of concentration in the 
compartment. Allowing vesicles to dwell between com- 
partments for a finite time causes the total number of 
molecules in the compartments to decrease, and hence 
yields an effective total concentration C c ff lower than the 
actual total concentration in the system C to t- 



N v 

Cos — Ctot — c* 



(10) 



Ja(C ' a )S % a (C ' a )P a ->p(C ' a , C p) 



(9) 



where N v is the number of vesicles between compart- 
ments, and CI the concentration carried by the i-th. vesi- 
cle. 

This sum over all the vesicles is actually a random 
variable, but its mean can be computed analytically in 
certain cases. To know where the symmetric/asymmet- 
ric transition occurs, we may compute C c g in a sym- 
metric system, and the result will be valid up until the 
transition. If the system is symmetric, vesicles have the 
same average concentration C v and Yli^i^v ~ N V C V . 
Let us assume each vesicle in the inter-cisternal medium 
has a rate of fusion W r towards any of the compart- 
ments. The average number of vesicles in the media is 
then 2K V /W r (where K v is the rate of individual vesi- 
cle secretion). Moreover, the average vesicle concen- 
tration C v can be written as the maximum concentra- 
tion C™ ax a vesicle may carry, times the average vesicle 
saturation fraction S (obtained from Eqj3|, leading to 
C cS = C tot - 2^C™ ax S. Finally, the product K v C™ ax 
is the number of vesicle leaving a compartment per unit 
time multiplied by the maximum concentration of each 
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vesicle, and can be identified with J s = KqC s . The crit- 
ical point of a system of total concentration Ctot with 
vesicles staying a finite time between the compartments 
can thus be obtained from the critical point (Eq[5]) of a 
system with infinitely fast fusion, but with an effective 
total concentration C^j™ given by: 



K 

^-fsym _ JV rr 

W r 



S = 



C ,, 



^syn 
°cff 



2C S 



(11) 



namely 



°eff 



C u 




with <j> s = Cg/C to t, w r = W r /K . 

Similarly, one can compute the effective concentration 
in a fully asymmetric system, in which one compartment 
has a concentration C^ yin and the other has a concen- 
tration close to zero: 



f-'eff - °tot 



W r 



S' = 



y^asym 
°eff 



^tasym 
°eff 



(13) 



The difference between Eq[l2]and Eq[l3]being that in the 
latter case, the empty compartment is sending out empty 
vesicles, and only vesicles from the first compartment 
contribute to the depletion effect. We find : 



£,asym | 



,W r ) = CT(l<Ps,W r ) (14) 



These predictions can be tested by numerical simu- 
lations (described in the Appendix). Comparison with 
the solution of the infinitely fast fusion model are shown 
in Figj2] Not only the location of the critical line, but 
also the actual values of the concentrations in each com- 
partment in the asymmetric steady state, where found to 
agree very well, even for low vesicle fusion rate. 



II. COMPARTMENT DIFFERENTIATION IN 
AN OPEN SYSTEM 

The relative simplicity of the model presented in 
Sec jl A| essentially characterized by two parameters 
(C s /C to t and C//C tot , Figjlja), allows us to address is- 
sues of direct biological relevance, such as the presence of 
external fluxes of material, and the possibility for chemi- 
cal transformation within the system. Organelles such as 
the Golgi are strongly polarized, with distinct entry and 
exit faces. We investigate the consequences of coopera- 
tive transport in such open systems assuming that a par- 
ticular species enters the system through compartment 
1, and exits through compartment 2, while exchange be- 
tween the two compartments proceeds as described pre- 
viously. Mathematically, this amounts to including a 
source term I± = Ji n and a sink term I 2 = —J out in 
Eq[T] yielding: 
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FIG. 2. Root mean square (RMS) of the difference of con- 
centration between the two compartments as a function of 
4>f — Cf/Ctot for various values of 4> s = C s /Ctot- Dash- 
dotted lines represent the mean-field values of the concen- 
tration (normalized by C^g ym /Ctot with effective parameters 
4>t f} = C 3 /C% m and <j> f = C//C^ m . Solid lines represent 
the simulated results with vesicles in the inter-compartments 
medium, with a vesicle fusion rate W r = Kq, i.e. up to 40% 
of the molecules are out of the compartments. The non-zero 
value of the RMS in the symmetric state is due to fluctuations. 



where a simple linear relationship was assumed for the 
out-flux: J out = K ?[C2, and where the exchange fluxes 
(^1^2) are still given by Eqs 2|3|4 At steady state, all 
fluxes must be balanced, including the net flux J between 
the two compartments: Ji n = J ou t = J. 



J=J X - 



J2 



d t C x = d t C 2 + K oS C 2 (15) 



A. Qualitative analysis 

The dynamical behaviour of the set of equations Eq[l5| 
is discussed in some details in Sec |IIB[ but a qualita- 
tive understanding of the open system may be inferred 
from the results obtained for a closed system. We showed 
(Fig|T|b) that the flux Ji^2 cannot exceed a maximum 
value and decreases upon increasing total concentration 
beyond a threshold. In an open system, this behaviour 
may result in the absence of a steady state: if the in- 
flux into compartment 1 exceeds the maximum net flux 
from 1 — > 2, the concentration C\ of the entry compart- 
ment steadily increases with time, leading to a further 
decrease of inter-compartment exchange. In the absence 
of other compensatory mechanisms, C\ would diverge 
and C2 would vanish, leading to a vanishing exit flux. 
This divergence is probably not realistic, but it illus- 
trates the consequence of such non-linear transport for 
an open system: beyond a critical influx, the system is 
essentially blocked, filtering transit proteins at a very low 
flux. While such feature has a negative impact on the rate 
of transport, it strongly increases the residency time of 
molecules and may prove advantageous to a system such 
as the Golgi apparatus, whose function is to process and 
chemically modify proteins. 
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FIG. 3. Phase-space trajectories of system with an exit flux J out = K aC2 (K a g = 0.005A"o), and an input flux Ji n = 0.005-K'oC'a 
(left) and Ji n = 0.025KoC s (right). Dash-dotted lines represent C2 = and dashed lines Ctot = 0. Red arrows represent initial 
condition with convergent trajectories whereas blue arrows are for initial conditions yielding a divergence of Ci. 



B. Phase-space trajectories of an open system 

We now discuss possible dynamical behaviours of an 
open systems satisfying the kinetic equation Eq |15[ where 
the fluxes between the two compartments are given by 
Eqs 2|3|4 Although the exchange rules between the com- 
partments are symmetric, the existence of external fluxes 
breaks the symmetry of the system, and different con- 
centrations should be expected in the two compartments 
even for low incoming flux. The critical behaviour at high 
incoming flux, as depicted in Fig[l] for a closed system, 
has nevertheless a profound impact on the steady states, 
or the absence thereof. 

As discussed above, one expects the flux exchanged 
between the two compartments to present a maximum 
value Jmax (necessarily smaller than the maximum pos- 
sible flux KqC s , see Figfljb), theoretically leading to a 
diverging concentration in the first compartment and a 
vanishing exchange flux if J in > J max - Depending on 
initial conditions, a non-convergent behaviour might ac- 
tually also appear for values of Ji n a priori compatible 
with the existence of a steady-state. For instance, if the 
initial concentration is very high in the first compart- 
ment, the divergent regime may occur for smaller in-flux 
Jin < Jmax ■ This can be understood by considering the 
phase space trajectories of the system. 

The coordinates in phase space are the concentrations 
(Ci, C2), and the steady states (if any) are given by the 
intersections of the Ci = and C tot = curves. Since 
>/out = K Q fiC2, the line Ctot = is obviously the line 
C2 = Jin/K ff, whereas the curve Ci = has to be 
computed numerically. If these two lines do not intersect, 
there is no steady state and Cj always diverges. If they 
do intersect, the thus-defined fixed points may be linearly 
unstable, or may be surrounded by a basin of attraction, 
as shown in Fig(3| 

The phase space representation show Fig|3]can be used 
to study the consequences of a transient change of the 
input flux (i.e. a pulse or a block of secretion). Let 
us consider a system which is in a stable steady state 
{C\,Ci) for an input flux J\ n . If the incoming flux is 



changed to Jf n at time t\, the phase space trajectories 
will be changed, and the system will follow a new tra- 
jectory starting from {C\,Cl). According to this new 
trajectory, the system will reach a new position (Cf, Cf ) 
at a time £2- If the flux is then switched back to its 
original value J} n , (Cf , Cf) will not necessarily be in the 
attractive region of the stable steady state. Therefore, 
a transient change of the incoming flux may bring the 
system out of a stable steady state. In the case of a 
strong pulse (Jf n » J} n ) the system may follow a di- 
vergent trajectory and the concentration C\ will increase 
strongly with time. Formally, whatever the (finite) value 
of (Ci,C|) after a pulse, the system may reach a sta- 
tionary regime if the incoming flux Ji n after the pulse is 
small enough. However, this may take a very long time. 
The approximation C\ — > 00 , J in = shows this time 
grows like {Ci) 2 . 



III. CONSEQUENCE OF COOPERATIVE 
TRANSPORT FOR PROTEIN MATURATION 
AND SORTING 

We now quantify the consequences of the kind of co- 
operative transport considered here on protein matura- 
tion and sorting. We investigate the situation sketched 
in Fig|4j where a molecular species A enters the system 
via compartment 1 and is transformed into a species B 
by maturation enzymes, before leaving the system via 
compartment 2. The processing accuracy is defined as 
the total fraction of the input that exits the system as 
mature (B) molecules: 



Accuracy : 



1 



(16) 



where A = J °° Ji n dt is the total amount of A molecules 
to have entered the system and J^ ut is the out-flux of B 
molecules. The accuracy thus defined reaches unity when 
no molecules exit the system without being processed 

(Jout = 0). 
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amount of time (which we assume very small), and set 
A\{t = 0) = A and Jj n = below. 

In order to focus on the role of cooperative transport, 
we further assume that particle export is not a rate- 
limiting step (K o g/K — > oo), and we analyze the pro- 
cessing accuracy in terms of a competition between the 
kinetics of maturation and transport (controlled by 4 pa- 
rameters) . 



FIG. 4. Sketch of an open system with protein maturation. 
Particle enter the system through compartment 1, undergo 
maturation A — > B while in the system, are exchange between 
compartment via cooperative transport, and exit the system 
via compartment 2. 



A Michaelis-Menten maturation kinetics is chosen in 
order to account for the limited amount of enzymes in 
the system. Calling A\ and B\ the concentrations of A 
and B in the first compartment, we have: 



d t B 1 = R{A\)A\ = R C V 



Ax 



Ai + G n 



(17) 



with an identical kinetics in compartment 2. Here, Rq 
is the maximal maturation rate and C m is the concen- 
tration of A beyond which enzymatic reaction saturates. 
For simplicity, we assume that the state (A or B) of a 
molecule influences neither its transport between com- 
partments nor its export from the system, so that Eq |15| 
is still valid for the concentrations C\, 2 = ^.1,2 + B\ 2 . 
Taking the weights of A and B in the fluxes to be their 
respective weights in the compartments : 



J A = 



J 1 



Ai 



Ai+B x 
Si 



A- 2 



Jl- 



A2+B2 

B 2 



J: 



2-»l 



Ai+Bi A 2 + B 2 

the following set of kinetic equations is obtained: 



A\ — Ji n 

A 2 = -R{A 2 )A 2 
B 2 = R{A 2 )A 2 



R{A 1 )A 1 
R(A 1 )A 1 



J B 



K oS A 2 



J u - K oS B 2 



(18a) 
(18b) 



(19a) 
(19b) 
(19c) 
(19d) 



Normalizing rates with the vesiculation rate Kq and 
concentrations with the concentration C s at which se- 
cretion saturates, Eq|l9] is controlled by 5 parameters. 
These are: ro = Rq/Kq and C m /C s , which compare the 
activity of the maturation enzymes and of the secretion 
machinery, Ct/C s , which defines the threshold for domi- 
nant specific fusion (Eqj5|, and K s/Kq, which compares 
exit and exchange rates. The fifth parameter is the nor- 
malized amount of material going through the system: 
Aq/C s . For simplicity, we investigate a situation similar 
to the so-called pulse-chase procedure |15j. where a fixed 
amount of material is delivered to the system in a finite 



A. Processing accuracy for linear transport 

In order to quantify the consequences of cooperativ- 
ity on the processing accuracy of a two-compartments 
system, we compute the accuracy of a perfectly linear 
system by linearizing Eqs 2|3|4 when A -C C m , C s ,Cf, 
yielding: J\^ 2 — K Ci/2 and J 2-^1 — K C 2 /2. Choos- 
ing a linear exit flux for simplicity (J ou t = K gC 2 , for 
which wc later take the limit K g — > 00), and the initial 
conditions d{t = 0) = 0.(0) and C 2 (t = 0) = 0, the ki- 
netic evolution of the vector C = {C\(t). C 2 (t)} is easily 
obtained: 



C(t) = e 



M,l 



O(o) 




Mt = — 



1 

1 1 



1 

2k oS 



(20) 



where fc g = _ ^ t - The matrix M; can be diagonalized, 
and the matrix exponential becomes a regular exponen- 
tial, and the concentration in the second compartment 
reads: 



C 2 (t) = 



0(0) 



2^1 



ft off 



with the eigenvalues 



a± 



K 



±Jl + k 2 oS -(l + k oS ) 



(21) 



(22) 



The (normalized) probability density that a particle 
exits the system from the second compartment at time t 
is Pexit{t)= M(t)/Ci(0): 



t(t) = 



Kok ff 



k 2 



') 



(23) 



The mean residence time of a particle in the system is 
thus (T) = / °° dt{tP exlt {t)) = 2(1/K + l/K oB ). 

The accuracy of protein maturation (A — > B) and sort- 
ing is defined as the fraction of the total quantity of 
molecules that entered the system to leave as matured 
molecules (Eq|l6|. It may also be written as: 



Accuracy = 



+00 



P exit (t)P(B,t\A,0)dt (24) 



where which P(B, t\A, 0) is the probability for a molecule 
to be mature (state B) at time t while starting imma- 
ture (state A) at t = 0. At the linear level, the mat- 
uration kinetics (Eq[l7| becomes: dtB a = RoA a , and 
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P(B,t\A,0) = 1 - e - Rot . The efficiency of the lin- 
ear system may then be computed analytically using 
Eqs.( |22|23|24| ), yielding : 



Accuracy | 



linear 



2r (1 + r + fcpg) 
k Q s + 2r (1 + r + k oS ) 



(25) 



with ro = Rq/Kq. Taking the limit fc ff — > oo as dis- 
cussed, the benchmark to which more complex transport 
and maturation processes must be compared is thus the 
linear accuracy: Accuracy |ii near — > 2r /(l + 2r ). 



B. Processing accuracy without specific vesicular 
fusion 



Saturation of maturation enzymes and transport in- 
termediates (C m ,C s < A <C C/, with A the initial 
particle concentration) has mixed effects on the systems 
processing accuracy. Saturation of inter-compartment 
transport at high concentration (for Aq C s ) causes the 
particle residency time of molecules to grow as A (Eq{3| 
while saturation of enzymatic reaction (for (A 3> Cm) 
causes the mean maturation time increase linearly with 
Aq (Eq[i~7|), so the net effect on processing accuracy de- 
pends on the precise values of the parameters. 

In order to get a fell for the role of the different param- 
eters, we compute the first order correction to the linear 
processing kinetics studied in Sec III A in the limit of 
very fast exit from the second compartment: AiofT ^ co- 
in this case, the accuracy is controlled by the flux exiting 
the first compartment, now written J out = \KqC s S(C\) 
and Eq[l9] may be rewritten: 



A = 

B = 



-C s 

-a. 



1 A 



C = A + B 
A 



2C + C : 
1 B 



roCr, 



A -\- C 71 
A 



2C + C S ^ r ° Crn A + C, 



(26a) 
(26b) 

(26c) 



where the subscript 1 has been dropped in the con- 
centrations, time has been normalized by 1/Kq, and 
r = Rq/Kq. Taylor expansion of this set of equation 
for Aq <C min(C s ,C m ) yields the first order correction 
to the accuracy of the linear system (Eq|25|: 



Accuracy = 



2r , A a C m (l + 2r )-C s (l + r ) 



1 + 2r C s 



O 



C m (l + r )(l + 2r f 

c s 



(27) 



An increase in processing accuracy at high concentra- 
tion requires that maturation saturates after secretion 
according to C m /C s > (1 + r )/(l + 2r ). This can be 
seen in Fig(5] which shows the variation of the system's 
processing accuracy as a function of the total amount of 
material to be processed, in the absence of cooperative 
fusion (dashed lines). 
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FIG. 5. Accuracy (Eq 16 1 as a function of the initial con- 
centration Aq, for different saturation ratios C m /C s of mat- 
uration and transport. At high concentration, specific vesicle 
fusion greatly enhance processing accuracy (solid lines, with 
Cf/Cs = 0.1), as compared to random fusion (C/ — > oo, 
dashed lines) (with Rq/Kq = 2, and K o «/K = 100). 



C. Processing accuracy and cooperative transport 

Cooperative fusion, when combined with saturation of 
the transport, leads to a robust increase of the processing 
accuracy of a compartmentalized organelle responsible 
for protein maturation and sorting (see Fig[5j solid lines). 
This increase can be understood as follows; at high con- 
centration (j4o > Cf,C s ), specific interactions promote 
backward fusion of vesicles secreted by the highly concen- 
trated compartment. As the forward fusion probability 
is very low (Pi-j.2 ~ EqQ the mean residency time 

increases as Aq, as compared to the linear incre ased ob- 
served in the absence of specific fusion (Sec III B I . On the 
other hand, the mean maturation time is still linear in 



Aq, so high concentrations lead to a more pronounced in- 
crease of the residency time compared to the maturation 
time, resulting in an increased processing accuracy at 
high concentration, even if the chemical transformation 
is performed by a limited amount of maturation enzymes 

(c m «a). 



IV. CONCLUSION AND OUTLOOK 

The predicted high processing accuracy displayed Figj5] 
essentially stems from the increase of the residency time 
of molecules transiting through the system. In strik- 
ing contrast with the usual Fick's law of gradient-driven 
transport, cooperative transport through the compart- 
mentalized system described here is strongly impaired by 
large concentration gradients. A strong prediction of our 
model is that the transport time actually increases with 
an increasing incoming flux (above a threshold). Pulse- 
chase experiments on the Golgi seem to show this trend, 
but data are still too scarce for a direct comparison (see 



Fig.4.Z in [H]). Although an apparent functional draw- 
back, slow transport through organelles is common. For 
instance, the typical transport time across the Golgi is of 
order of 20 minutes [IS] , whereas diffusion of a membrane 
protein over an area equal to that of the entire Golgi ap- 
paratus (of order 10/im 2 ) should be of order one minute 
(with a diffusion coefficient ~ 0.1/xm 2 /s [I]). 

In this paper, we showed that organelles constantly 
exchanging material via transport vesicles may sponta- 
neously adopt different biochemical identities, provided: 
i) the flux of vesicles secreted by an organelle is bounded, 
and ii) there exists some level of specific vesicle-organelle 
fusion directed by molecular recognition. In open sys- 
tems traversed by fluxes, these transport properties give 
rise to a dynamical switch from a linear to a low through- 
put kinetics above a critical influx. For compartmental- 
ized organelles whose function is to process and export 
influxes of proteins, such as the Golgi apparatus, this 
switch allows the export rate to spontaneously adjust to 



the amount of material to be processed, a definitive func- 
tional advantage that may avoid the release of unpro- 
cessed material even under high influx. Future extension 
of the present model to multi-component transport (see 
S.I.) will allow to assess the importance of the exchange 
of membrane area (lipids) between compartments (the 
compartment sizes were assumed constant here). It will 
also allow us to explore the role of specialized membrane 
domains; domains for protein processing and domains for 
protein export, which have recently been reported |16) . 
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Appendix A: Numerical simulation - Finite vesicle fusion time 

To test the predictions of Eqs ]13|14| in the case vesicle fusion occurs with a finite rate, we performed a numerical simulation 
of a system with finite vesicle fusion time and a total concentration Ctot and compared the location of the critical line with 
the infinitely fast fusion model, the equations of which we solved numerically. The numerical simulation consists of two 
compartments of concentrations Ci and C2 from which vesicles may bud with a rate K v . Each vesicle budding from a 
compartment a has a saturation S(C a ). At each timestep, each vesicle may merge with a compartment at a rate W r , and the 
compartment is chosen according to the probability P described in the main text. The algorithm may be written as follows : 

#Fusion prob ability of a vesicle of saturation Sv with the first compartment 
def Pfl(Sv,Cl,C2) = (Sv*Cl+Cf)/( 2*Cf + Sv*(Cl+C2) ) 

#Saturation of the vesicles leaving from a compartment of concentration C 
def S(C) = C / (C+Cs) 

#SvesfiJ: Saturation of the i—th vesicle 
#Nves : number of vesicles 

while t<Tmax : 
t=t+dt 

#Checks if a vesicle leaves the first compartment 
if rand (1) < Kv*dt : 
Nves=Nves+l 
Sves [Nves] = S(Cl) 
C1=C1 - Cv*Sves [Nves] 

#Checks if a vesicle leaves the second compartment 
if rand (1) < Kv*dt : 
Nves=Nves+l 
Sves [Nves] = S(C2) 
C2=C2 - Cv* Sves [Nves] 

#Checks for each vesicle if it merges with a compartment 
for i=l to Nves : 

if rand (1) < Wr*dt : 

if rand(l) < Pfl ( Sves [ i ] , CI , C2) : 
Cl=Cl+Cv*Sves [ i j 

else : 

C2=C2+Cv*Svcs [ i ] 
Sves [i]=0 
reorder (Sves , Nves ) 



